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Using density-functional theory we assess the stability of bulk and surface oxides of the late Ad 
transition metals in a "constrained equilibrium" with a gas phase formed of O2 and CO. While the 
stability range of the most stable bulk oxide extends for ruthenium well into gas phase conditions 
representative of technological CO oxidation catalysis, this is progressively less so for the 4d metals 
to its right in the periodic system. Surface oxides could nevertheless still be stable under such 
conditions. These thermodynamic considerations are discussed in the light of recent experiments, 
emphasizing the role of (surface) oxides as the active phase of model catalysts formed from these 
metals. 
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I. INTRODUCTION 

Rhodium and palladium exhibit a high reactivity 
for the CO oxidation reaction, which has led to their 
widespread use in catalytic car exhaust convertersi. The 
degree of oxidation of these materials in the reactive 
environment has been a frequent object of speculation 
i n th e e xpe rimental catalysis literature (see e.g. refs. 
I2l l3l4l l5ll6l l7l8l) . with oxide formation at the surface of the 
metal particles potentially having either beneficial or 
detrimental effects on the overall activity. Lacking mi- 
croscopic in-situ data to really resolve this issue a pos- 
sibly beneficial role has e.g. been suspected for PdO at 
Pd7, whereas Rhodium oxides have predominantly been 
viewed as poisoning the catalytic reactionS*^. 

While the atomic-scale techniques of Ultra High Vac- 
uum (UHV) surface-science would in principle be able 
to provide additional insight into this subject matter, 
formation of oxides in particular of the more noble met- 
als was often neglected in corresponding works, partly 
due to the fact that these oxides form only at rather 
high oxygen pressures and elevated temperatures. This 
emphasis on metallic substrates has changed over the 
last years, bringing oxide surfaces now also on the 
agenda of surface-science studies attempting to bridge 
the pressure gap between UHV and technological ox- 
idation catalysis^. Concomitantly, on Ru(0001)i2*ii, 
Pd(100)^^>i^, and Ag(lll)i^ii^ oxygen-rich environmen- 
tal conditions were reported to lead to oxide (or "surface 
oxide" ) formation, in all cases connected with a signifi- 
cant increase in catalytic activity. 

This new focus on the effect of the environment on 
surface structure and reactivity is intertwined with cor- 
responding theoretical efforts. In particular the concept 
of first-principles atomistic thermodynamica^^'i^'i^'^ has 
recently proven to be most valuable for the study of oxide 
surfaces (see e.g. refs. l2(Tll2l2l2l2l . In this branch of 



thermodynamics density-functional theory (DFT) is em- 
ployed to compute free energies, in order to identify the 
lowest-energy atomic structure for a given condition of 
thermodynamic reservoirs representing the surrounding 
gas phase, and possibly combining this with constraints 
to facilitate the analysis22iS4. In the following we will 
use and describe this formalism to discuss the stability 
of bulk and surface oxides of the late id transition met- 
als (TMs) Ru, Rh, Pd, and Ag. The results of these very 
basic thermodynamic considerations are found to nicely 
embrace the aforedescribed existing experimental data, 
e.g. pointing at the relevance of thin surface oxide layers 
on Pd under the conditions of high-pressure CO oxida- 
tion catalysis. 



II. "CONSTRAINED THERMODYNAMICS" 

The question we like to ask in this context is, given 
a certain gas phase environment composed of O2 and 
CO, does a bulk or surface oxide of the late 4d TMs 
form a stable phase? In a thermodynamic description we 
then resort to some form of equilibrium that exists be- 
tween various components i in the system, all of which 
are present in sufficient quantities. Under conditions of 
constant temperature, T, and pressures, {pi}, this al- 
lows to discuss relations between the respective reser- 
voirs, characterized by their chemical potentials fj,i(T^pi) 
or Gibbs free energies gi (T, pi) . For a gas phase composed 
of CO and O2 a forthright assumption of full equilibrium 
does, however, not yield a description suitable for our 
interest in oxidation catalysis: From energy considera- 
tions alone, CO2 would then result as the most stable 
gas phase molecule for almost all temperature and pres- 
sure conditions. On the other hand, because of the large 
free energy barrier for the gas phase reaction CO -I- 1/2 
O2 — > CO2, such a full equilibrium will not be attained on 
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time scales relevant to us. Instead, we may ignore CO2 
formation in the gas phase and consider the presence of 
two separate, independent reservoirs for O2 and CO in 
the environment. While both components are therefore 
not in equilibrium with each other, each one is individu- 
ally assumed in equilibrium with either the metal (M) or 
the oxide (MxOy) phase. 

As a consequence the chemical potentials of both oxy- 
gen and CO in the whole system are then determined by 
the surrounding gas phase reservoirs, i.e. their tempera- 
ture and pressure dependence is given b y^^'^^ 



/^o(7',P02) = ^ 



and 



total 



Ao.(T,/) + fcBm^ 



(1) 



^^Co{T,pco) = E'cT' + Mco(r,/) + knTln (^^^ , 

(2) 

where the temperature dependence of p.o^iTjp'^) and 
MCo(r,p'^) includes the contributions from vibrations 
and rotations of the molecules, as well as the ideal gas 
entropy at p° = 1 atmosphere. The latter quantities can 
be computed from first principles, yielding results that 
are at the temperature range of interest to us virtually 
indistinguishable from the experimental values listed in 
thermodynamic tables^^iS^. Via eqs. and l(2Jl our re- 
sults obtained as a function of the chemical potentials, 
i.e. in (/io, Mco)-space, can then be converted into pres- 
sure scales at any specific temperature. Roughly repre- 
senting the temperature range relevant to our study, this 
will be illustrated with specific scales at T = 300 K and 
T = 600 K in all figures below. 

In the resulting "constrained equilibrium" 
situation^i^ the only pathway to CO2 formation 
is due to a reduction of the oxide or due to the catalytic 
CO oxidation at the substrate surface. In a flow reactor 
this constant build-up of CO2 is counteracted by the 
continuous stream (and thus removal) of gases over the 
catalyst surface, yielding a CO2 concentration gradient 
under stationary operation conditions. Although formed 
CO2 molecules may in principle readsorb at the surface 
(dissociatively or as a molecule), the probability of this 
event is very low compared to the much more frequent 
O2 and CO (re)adsorption. The surface is therefore 
unlikely to equilibrate with the surrounding CO2 gas 
phase, and a proper description would have to involve 
a dynamical treatment of the CO2 gas flow close to 
the catalyst surface. Such treatment will be published 
elsewhere2&. Here we use a very simple estimate for 
the CO2 free energy, namely the internal energy of 
a free CO2 molecule. The error introduced by this 
estimate will vary from material to material, but in 
our preliminary kinetic Monte Carlo simulations for 
RuO2(110) we find it to be small»2£ 

Although crude, this approximation is appropriate for 
the present trend study aiming only at a first assessment 
of oxide stability, and an uncertainty of some tenths of eV 



does not change the nature of our conclusions. In this re- 
spect we stress that the whole concept of a "constrained 
equilibrium" is not designed to yield precise answers to 
the microscopic state and reactivity of the catalyst sur- 
face. For that, the very dynamic behavior must be mod- 
eled by statistical mechanics. Yet, with atomistic ther- 
modynamics one can trace out the large scale behavior 
and identify those conditions where such a more refined 
treatment is necessary. 



III. BULK OXIDE STABILITY 

If there is a pure oxygen environment, the condition 
for the stability of the bulk oxide is 
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Here A^o is defined as ^io - (l/2)£'^°*'*'. For T = 
K the bracket on the right hand side equals half of the 
low temperature limit of the heat of formation, Hf{T — 
OK,p ~ 0). Even for higher temperatures this will 
not change substantially as for bulk phases the {T,p)- 
dependence of their Gibbs free energies is rather small, 
i.e. the right hand side will still be rather close to the 
Hf{T — QK,p = 0)-value and not follow the pronounced 
variation of Hf (T, p) (which is primarily due to the (T, p)- 
dependence of the /xoa-term, while only Eq^^^ enters the 
right hand side of eq. Q). We therefore replace the 
bracket in eq. Q by (l/y)iJ/(T = OK) and arrive at 
our first stability condition 



A^o ~ -Hf{T = OK). 

y 



(5) 



Additionally, the oxide can also be destroyed (reduced) 
by carbon monoxide. In a pure CO environment, the 
stability condition for the oxide is 



5m" o„ + y A*co 



< 
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y A^COa 
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Following the discussion above we will crudely approxi- 
mate /ico2 ill the present study by the internal energy of 
a free CO2 molecule. In a similar fashion as for the pure 
oxygen environment one can then simplify the stability 
condition to 



A^co^---ff/(T = OK) 

y 



AE' 



mol . 



(7) 



E^Q-^ and using the short 



introducing A^co = ^^co 

hand A£;™°i- = (E^^St " ^co'^ - 1 /'^E}^'^'^) for the differ- 
ence in binding energies of the three gas phase molecules. 
As described in detail in a previous publication^'* our 
DFT computations including zero-point vibrations give 
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TABLE I: Low-temperature limit of the heat of formation, 
AHf{T = OK), of the most stable bulk oxide of the four late 
4d transition metals from DFT and from experiment. The 
DFT values do not include the zero-point energies of the solid 
phases. 





Theory 


Experiment 


Ru02 


-3.4eVi22i 


-3.19e\42i 


RhzOs 


-3.8eVi2a 


-3.57e\42i 


PdO 


-O.geViS 


-0.88 eVi2i 


AgaO 


-0.2 eVi^ 


-0.34 eViSi 



^^moi. _ — 3 2eV, which agrees well with the exper- 
imental value of — 2.93eV~. Although state-of-the-art 
DFT total energies of the O2 molecule are subject to 
a noticeable error, some cancelation apparently occurs 
for the composite quantity Ai?™°'-, rendering it rather 
unimportant on the scale of the present study whether 
the experimental or theoretical value for AE"^"^- is used. 

Finally, if O2 and CO are both present in the gas phase 
in "constrained equilibrium" with the oxide, then the 
general stability condition is simply obtained by com- 
bining eqs. (O and Q for the separate cases: 

A/ico - A^o ^ = OK) + Ai?'"°S (8) 

while eq. (jSJ still applies for very low CO concentrations. 
Hence, in both conditions for the stability of the bulk ox- 
ide, eqs. ((3) and 0, the only quantity left to determine 
is AHf{T = OK). This value has been computed in pre- 
vious worksi^iSSiSSi^s and we content ourselves here with 
simply listing the DFT values for the most stable ox- 
ides of the four late 4d TMs in Tabled again comparing 
with the low-temperature limit of the experimental heat 
of formationSi. Keeping the sizable error in the O2 DFT 
total energy in mind, we notice in all cases a fortuitous er- 
ror cancelation, leading again to rather small differences 
when evaluating the stability conditions with either ex- 
perimental or theoretical numbers. Only in the case of 
silver oxide has AHf{T = OK) become so small that this 
error becomes somewhat disturbing (in particular for the 
stability condition in eq. ©)■ More refined treatments 
for approximating the stability conditions, e.g. including 
entropic contributions to AHf{T,p), can be devisedi^, 
but for our present purposes it suffices to see that bulk 
Ag2 has such a low stability that it will not play a role 
in our further discussion. 

The derived (T, p)-ranges for stability of the 4d bulk 
oxides are shown in Fig. ^ exhibiting a pronounced de- 
crease from Ru02 towards the oxides of the metals to its 
right in the periodic system. We stress that the aforedis- 
cussed crude treatment of the CO2 free energy translates 
into rather large error bars on the diagonal lines running 
from bottom left to top right in the figure (i.e. the limit 
due to the stability condition of eq. lO). The uncertainty 
in /ICO2 translates directly into the height of these lines, 
which with the present error estimate could therefore be 
shifted up or down by some tenths of an eV. 
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FIG. 1: Stability regions of bulk oxides of the late 4d transi- 
tion metals in (A/io, A/^co)-space, on the basis of the DFT 
data listed in Table L Additionally, pressure scales are drawn 
for T = 300 K and T = 600 K. The white line indicates con- 
ditions corresponding to po = Pco = 1 bar and 300 K< T < 
600 K, and environments relevant for CO oxidation catalysis 
over these metals would correspond to the near vicinity of this 
line. More specifically, the small white area indicates roughly 
the (r, p)-conditions employed in a recent experimental study 
by Hendriksen and Frenken (see text)— i— . The stability of 
Ag2 is so low, that is does not appear any more within the 
limits considered here. We stress that the rough treatment 
of the CO2 free energy translates into an uncertainty in the 
position of the diagonal lines (which could be shifted up or 
down by some tenths of an eV). 

IV. SURFACE OXIDE STABILITY 

Defering the discussion of the derived bulk oxide sta- 
bility to the end of this paper, we proceed to analyze in 
the stability context a second, only recently emphasized 
aspect of metal oxidation, namely the formation of thin 
surface oxides. While traditionally such films were con- 
ceived as closely-related thin versions of the correspond- 
ing bulk oxides, recent atomic scale characterizations of 
initial few-atom thick oxide overlayers especially on Pd 
and Ag surfaces revealed structures that had only little 
resemblance to their bulk counterparts, and/or were to 
a large degree influenced by a strong coupling to the un- 
derlying metal substratei^i^fiiSl. Due to this coupling and 
structures particularly suited for layered configurations, 
one may expect the stability range for such surface oxides 
to exceed that of the hitherto discussed bulk oxides. 

To determine the range of (T, p)-conditions in which 
a surface oxide would represent the thermodynamically 
most stable state, we follow the approach of Li, Stampfl 
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and Scheffler— , i.e. we evaluate the Gibbs free energy of 
adsorption and compare it to other possible states of the 
system like on- or sub-surface oxygen phases or the bulk 
oxide. For the case of a pure O2 environment and using 
the same approximations as discussed in the last section, 
we may write this Gibbs free energy of adsorption as 

AG(Aa*o) « 
« -\ (^^oAm - - No l^-E'C' + Amo)) 

- ^ {E^S^M + ^t^o) ■ (9) 

Here, i^o@M ^'^'^ ^m*'*' are the total energies of the sur- 
face with and without oxygen coverage of A^^o oxygen 
atoms per surface area A. In the second line, we have 
identified the first terms in the brackets with the aver- 
age binding energy of oxygen in the particular surface 
configuration and with respect to 1/2 Eq'^'^^. 

Obviously, the higher the oxygen content of a consid- 
ered surface structure, the steeper its AG'(A^o) will de- 
crease with increasing chemical potential. In the limit- 
ing case of an infinitely thick bulk oxide on top of the 
metal substrate, this will result in a vertical line that 
crosses the zero-axis at the stability condition for the 
bulk oxide in eq. For any higher A/^o the bulk ox- 
ide will be the stable phase, but the interesting question 
is to see whether there are lower A/iQ than this limit, 
for which the lines of a surface oxide structure turn out 
lower than e.g. on-surface adsorption or the clean surface 
corresponding to AG(A^o) — 0. 

We exemplify this ansatz by first analyzing the stabil- 
ity of a surface oxide on Pd(lOO). In addition to the well 
known p{2 x 2) and c(2 x 2) ordered adlayers with O in 
fee sites^i^, an ensuing {s/b x \/b)R21°-0 surface ox- 
ide phase on this surface was recently characterized as a 
rumpled, but commensurate PdO(lOl) film with a strong 
coupling to the underlying Pd(lOO) substrata^ii^S. Eval- 
uating the reported binding energies in this study with 
eq. © we obtain the results displayed in Fig. |5| Inter- 
estingly, the lS.pLo-ia,Tige where this surface oxide exhibits 
the lowest Gibbs free energy of adsorption not only ex- 
ceeds the one of the aforediscussed PdO bulk oxide, but 
also the range where the on-surface adsorption phases 
are more stable than the clean Pd(lOO) surface. In other 
words, the latter on-surface phases never correspond to a 
thermodynamically stable phase, and their frequent ob- 
servation in UHV experiments^i^ appears to be a mere 
outcome of the limited O supply offered, as well as of 
kinetic barriers e.g. for O penetration at the low tem- 
peratures employed (UHV experiments are typically per- 
formed by depositing a finite number of adatoms, rather 
than by maintaining a given gas pressure). 

It is interesting to compare these findings with the 
equivalent situation on another Pd surface orientation. 
On Pd(lll), a surface oxide not resembling the bulk- 
structure of PdO at all was recently reported as the 
product of high oxygen exposure and following the well- 
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FIG. 2: Computed Gibbs free energy of adsorption for the 
p(2 x 2) and c(2 x 2) on-surface adsorption phase, as well 
as for the (^5 x ^)_R27°-0 surface oxide on Pd(lOO). The 
surface unit cell of Pd(lOO) is 7.8 A^. The stability ranee of 
the surface oxide extends well beyond that of the bulk oxide, 

given by A/xo ^ — 0.9eV, cf. eq. (|KJ. The dependence on 
A/io is again translated into pressure scales at T = 300 K 
and T — 600 K for clarity. In the bottom of the figure, the 
"material type" which is stable in the corresponding range 
of O chemical potential is listed and indicated by the shaded 
regions. The binding energies used to construct this graph 
are taken from refs. i31ii32l 



studied p{2 x 2) on-surface adsorption phasc'^". Using the 
DFT binding energies provided in this study to evaluate 
eq. © , we again obtain a range of (T, p)-conditions out- 
side of the bulk oxide stability range, for which this sur- 
face oxide phase is most stable as shown in Fig. This 
time, however, this range extends only within -1.2 eV 

^ A/iQ ~ — 0.9eV, below which first the on-surface ad- 
sorption phase and then the clean Pd(lll) surface be- 
comes more stable. If we compare this to the results 
shown in Fig. |21 there are thus environmental conditions, 
in which a surface oxide may already be thermodynami- 
cally stable on Pd(lOO), while on-surface adsorption still 
prevails on Pd(lll). Knowing that these two surfaces 
form the predominant surface area of Pd nanoparticlea^^, 
this may have profound consequences on the oxidation 
behavior of the latter and we will return to this point 
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FIG. 3: Computed Gibbs free energy of adsorption for the 
p(2 X 2) and hypothetical p(2 x 1) on-surface adsorption phase, 
as well as for the so-called (\/6 x \/6)-0 surface oxide on 
Pd(lll). The surface unit cell of Pd(lll) is 6.7 Again, 
the stability range of the surface oxide extends beyond that of 
the bulk oxide, cf. Fig. \^ but this time there is also a finite 
range in which an adlayer forms the thermodynamically most 
stable phase. The binding energies used to construct this 
graph are taken from ref. ISOl 



below. 

While a similar situation with the commensurable 
p(4 X 4) surface oxide^** slightly exceeding the stabil- 
ity range of bulk Ag20 was also found for Ag(lll)^^, 
no (T, p)-range most favorable for the various intermedi- 
ate precursors suggested in the oxidation pathway from 
Ru(OOOl) to RuO2(110) is obtained when evaluating the 
DFT binding energies reported in the corresponding refs. 
l3&.37> Tentatively, we therefore expect surface oxides to 
play a thermodynamic role more for the 4d TMs towards 
the right of the periodic system, where the decreasing 
thermal stability of the bulk oxides and the lowered bulk 
modulus of oxide and metal phase enhances the influence 
of oxide-metal coupling, and thus the tendency to form 
commensurable, even non-bulk like surface oxide config- 
urations. It is interesting to notice that the extended 
stability range of the surface oxide phase discussed on 
Pd(lOO), cf. Fig. 13 corresponds more or less to the 
range exhibited by bulk Rh203 in Fig. ^ i.e. it comes 
now rather close to the schematically indicated range 



of (T,po,Pco)-conditions representative of technological 
CO oxidation catalysis. 



V. IMPLICATIONS FOR THE CATALYTIC CO 
OXIDATION 

Due to the particular construction of our "constrained 
equilibrium" approach we could expect the results sum- 
marized in Figs, n to 121 to give a rough reflection of the 
state of the catalyst during steady-state CO oxidation, as 
long as kinetic effects both in the formation of the (sur- 
face) oxides themselves, as well as in the ongoing reac- 
tion on the surface (e.g. possibly consuming oxygen more 
rapidly than it is replenished) are neglig ible2&. Although 
technological catalysis typically operates well above room 
temperature, such effects can never be ruled out in gen- 
eral. Clearly, a comprehensive understanding of the full 
catalytic activity and the microscopic state of the cata- 
lyst (i.e. its very dynamic behavior) requires modeling by 
statistical mechanics, and the results of thermodynamic 
considerations must consequently be treated with great 
care. Still, we believe that our results may be employed 
for a first qualitative assessment of oxide formation on 
the late 4c? TMs, in particular if supplemented by some 
already known kinetic constraints. 

In this regard we first notice the rather large stability 
range of bulk Ru02, extending well over the whole range 
of environmental conditions relevant for technological CO 
oxidation catalysis, cf. Fig. ^ Rather high temperatures 
of the order of T ~ 600 K are required to form this bulk 
oxide^"'^"'^ , which would explain the recently reported lack 
of oxide formation at Ru(OOOl) even under high pressure 
conditions at temperatures around T ~ 400KiiS. Yet, 
once this bulk oxide is formed, our thermodynamic re- 
sults imply that the catalytic activity of Ru catalysts 
might in reality be due to oxides, in line with the results 
of a number of recent experimental studies^iiSiii. 

For Rh the stability range of its bulk oxide has al- 
ready decreased significantly. Even if the catalyst was 
in a situation close to thermodynamic equilibrium, our 
results suggest that either metal or bulk oxide could pre- 
vail at technologically relevant environmental conditions, 
depending on the exact partial pressures and tempera- 
ture in the gas phase, i.e. depending on whether one is 
below or above the stability line in Fig. ^ (and recall 
that the present uncertainty in the location of this line 
is some tenths of eV). Even more intriguing, oscillations 
between the two phases could result as a consequence of 
fluctuations, when operating under environmental con- 
ditions very close to the stability line of Rh203. Simi- 
lar considerations might also be possible for Pd, yet this 
time not concerning transitions to bulk PdO (which is 
already far too unstable), but to thin oxide films lim- 
ited to the very surface region - as mentioned in section 
IV the stability of the surface oxide on Pd(lOO) extends 
over a similar range to that of bulk Rh2 03 in Fig. ^ 
In this respect we have also roughly marked in Fig. ^ 
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the (T, p)-conditions employed in a recent experimental 
study on Pd(lOO), where a substantial roughening of the 
surface with presumed constant formation and reduction 
of ultra-thin oxide domains during the CO oxidation re- 
action was reported^^'^'^. The closeness of these experi- 
mental environmental conditions with our rough instabil- 
ity estimate for the surface oxide on Pd(lOO) is startling, 
in particular when further recalling the frequent observa- 
tion of spatio-temporal pattern formation during the CO 
reaction at other Pd surfaces like Pd(llO), which was 
there also connected to reversible oxide formation and 
reduction^. 

At this stage it is worth digressing on the fact that such 
thin oxide films with possibly even fluctuating conditions 
at the catalyst surface are tough to describe (even con- 
ceptually) within the traditional language of generic reac- 
tion mechanisms (see e.g. ref. [s^for a review). For metal 
surfaces mostly a Langmuir-Hinshelwood (LH) type be- 
havior is discussed, where both rcactants adsorb at the 
surface first and react thereafter. For the case of oxida- 
tion reactions on oxide surfaces on the other hand. Mars 
and van Krevelen (MvK)i2i suggested that the underlying 
substrate itself might become continually reduced and 
(re-)oxidized in the on-going reaction, i.e. that oxygen 
from the oxide lattice is consumed in the reaction and 
then replenished from the gas phase. Within the frame- 
work of macroscopic rate equations, where this approach 
was made and has ever since been very popular (in the 
oxide literature) , this translates to relaxing the constraint 
of a fixed number of active surface sites underlying the 
LH kinetics. The resulting MvK rate equation is then 
similar, but not identical to the one derived for a LH re- 
action mechanism, although in actual practice most data 
can be fitted equally well with both rate equationgSS. 

This only vague distinction between both mechanisms 
becomes even more blurred when starting to think mi- 
croscopically and particularly when considering what is 
now called a surface oxide (a 1-2 layer thick oxide-like 
film). For the case of a fixed oxide lattice, where cer- 
tain surface oxygen atoms just take part in the reac- 
tion and the resulting vacancies are refilled afterwards, 
a distinction between LH and MvK is in fact just seman- 
tics: One could equally well view the oxide lattice with 
the created vacancies (thus MvK) as the real substrate, 
where then oxygen adsorbs into these sites and reacts 
thereafter (thus LH). Although such a scenario has been 
advocated as a MvK mechanism for CO oxidation on 
RuO2(110)iS, one might even argue that it contradicts 
the original MvK definition in its strictest sense, which 
emphasized the non-finite number of active sites (that is 
certainly fixed by the number of surface oxygen sites in 
the latter situation). 

Other interpretations of MvK have therefore focused 
more on a diffusional aspect that would arise if oxygen is 
consumed at one point of the lattice surface but replen- 
ished elsewhere, requiring some form of oxygen transport 
in or at the surface (that might even become rate-limiting 
and thus yield a new type of kinetics)^^. Nevertheless, a 



clearcut distinction could even then not always be made. 
Picturing in particular a possibly only one-layer thick 
oxide film, that does not necessarily resemble a bulk-like 
oxide structure at all, or the aforedescribed situation of a 
surface close to an instability with a continuous formation 
and destruction of a thin surface oxide (possibly creating 
mesoscopic roughness due to the etching of metal parti- 
cles from the oxide framework^^'S) , how would one be 
able to distinguish what is an oxide vacancy, what the ox- 
ide, what an oxygen atom diffusing on or below the metal 
surface, let alone determine whether this would show up 
in a unique kinetics that would still be describable by 
simple macroscopic rate equations? 

In this regard, another intriguing result of our thermo- 
dynamic considerations is the different stability range of 
surface oxides on Pd(lll) and Pd(lOO). Again, provided 
the catalyst is close to thermodynamic equilibrium, this 
would imply that under certain (T, p)-conditions the ac- 
tive phase for CO oxidation could be on-surface adphases 
(i.e. LH) on some facets of nanoparticles and surface ox- 
ides (let's call this MvK) on others. A situation that 
could make the overall kinetic data (and its modelling) 
much more complex than hitherto assumed. In this par- 
ticular sense, an understanding of the data from Pd could 
be much more difficult than say for Ag, where we would 
not expect neither bulk nor surface oxides to play a sig- 
nificant role in CO oxidation catalysis on the basis of the 
present results, and in agreement with the conclusions of 
a more detailed recent theoretical studjsi^. 



VI. SUMMARY 

In conclusion we have derived simple thermodynamic 
stability conditions for bulk and surface oxides in contact 
with an environment formed of CO and O2. Applying 
these to the late Ad transition metal series from Ru to 
Ag, we obtain a progressively smaller stability range for 
the bulk oxides that extends well into the environmental 
conditions representative of technological CO oxidation 
catalysis only in the case of RUO2. In particular for Pd, 
thin surface oxides could nevertheless still be stable in the 
more oxygen-rich part of this (T, p)-range, while for CO 
containing gas phases neither surface nor bulk oxide will 
be stable on Ag. Most intriguingly, oscillations between 
metal and bulk (surface) oxide phase for Rh (Pd) are con- 
ceivable under technologically relevant (T,p)-conditions. 

Provided the catalyst is close to the constrained ther- 
modynamic equilibrium forming the basis of our consid- 
erations, these results may give first insight into the role 
of oxide formation played in technological CO oxidation 
catalysis at these 4d TMs. A comprehensive understand- 
ing of the microscopic state and full catalytic cycle must 
however be based on statistical mechanics as the next 
step, i.e. by modelling the very dynamic behavior and 
interplay of the manyfold of elementary processes. 
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